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Abstract. The temporal properties of a sample of 498 bright gamma-ray bursts (GRBs) with durations between 
0.05 s and 674 s were analysed. The large range in duration (Tgo) is accompanied by a similarly large range in 
the median values of the pulse timing properties including rise time, fall time, FWHM and separation between 
the pulses. Four timing diagrams relating these pulse properties to Tgo are presented and show the power law 
relationships between the median values of the 4 pulse timing properties and Tgo, but also that the power laws 
depend in a consistent manner on the number of pulses per GRB. The timing diagrams are caused by the correlated 
properties of the pulses in the burst and can be explained by a combination of factors including the Doppler boost 
factor r, a viewing effect caused by a jet and different progenitors. GRBs with similar values of Tgo have a wide 
range in the number of pulses. GRBs with the large number of short and spectrally hard pulses may occur either 
from a homogeneous jet with a higher average value of F or close to the axis of an inhomogeneous jet with higher 
values of F near the rotation axis. The less luminous GRBs with fewer pulses may originate further from the axis 
of the inhomogeneous jet. The pulses in GRBs have six distinctive statistical properties including correlations 
between time intervals, correlations between pulse amplitudes, an anticorrelation between pulse amplitudes and 
time intervals, and a link to intermittency in GRS 1915+105. The timing diagrams and correlated pulses suggest 
that GRBs are powered by accretion processes signalling jets from the formation of black holes. 
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1. Introduction 

The cosmological origin of gamma-ray bursts (GRBs) re- 
quires an extraordinary amount of energy to be emitted 
in gamma-rays ( Piran, 1999 ). The source of this energy 
is thought to be a cataclysmic event involving mergers of 
compact objects such as neutron stars or a neutron star 
and black hole ( Paczynski, 1991 ; Ruffcrt fc Janka, 1999 ) or 



the formation of a black hole in massive stars (MacFadyen 
fc Woosley, 1999| ; |Popham et al., 1999|). GRB light cu rves 
are complex and irregular ( Fishman fc Mccgan, 1995| ) and 



2. Pulse analysis in GRBs 

The BATSE experiment has provided the most uniform 
collection of GRBs available. A large sample of the bright- 
est GRBs with data combined from the four energy chan- 
nels was analysed using either wavelets or median filters. 
A detai led account of this pro c ess has been described else- 
where dQuilligan et al., 2002| ; |McBreen et al., 200 1| ). An 



a range of techniques have been developed to elucidate 



their structure 


McBreen et al., 1994; 


Morris et al., 1996; 


Lee et al., 2000; 


Beloborodov et al., 2000; 


Quilligan et al., 



2002). In this paper we present the results of the analy- 
sis of pulses in a large sample including long and short 
GRBs. The main results are presented in section 3 includ- 
ing four timing diagrams and our analysis of pulse prop- 
erties. These results along with the unique properties of 
the pulses in GRBs are interpreted in section 4 and used 
to provide additional insights to the emission process and 
the central engine that support jets from accretion into 
newly formed black holes. 
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automated pulse selection algorithm was used to detect 
pulses in the denoised signals. The samples included the 
brightest 319 GRBs with Tgo > 2 s analysed at 64 ms 
resolution and 100 GRBs with Tgo < 2 s analysed at 5 ms 
resolution. To extend the sample to include GRBs with 
Tgo well beyond 100 s, a further 79 bright GRBs were 
included and analysed at 256 ms resolution. 

Pulses that were isolated from neighbouring pulses at 
>50% level were selected and the rise times i r , fall times 
tf and FWHM were measured. The time intervals between 
pulses AT, above the 5<r threshold, were determined. The 
individual distributions of t r , tf, FWHM and AT have 
a wide range and are consistent with lognormal distri- 
butions. The bursts were split into Tgo categories with 
comparable number of pulses and typical lower bound- 
aries 0, 0.5, 2, 15, 30 and greater than 100 s for the GRBs 




Fig. 1. The timing diagrams in GRBs. The median values obtained from the lognormal distributions for the pulse 
timing parameters (a) t r , (b) tf, (c) FWHM and (d) AT are plotted versus duration T90 for GRBs in three categories 
i.e. 1 < N < 3 (blue), 4 < N < 12 (red) and N > 12 (green). The crosses signify the range covered in T90 and e /i±CT for 
the lognormal distribution which includes 16% to 84% of the pulse values for that T90 bin. The upper diagonal line is 
the limit where the pulse parameter is equal to T90 and the lower lines by the limited time resolution i.e. 5 ms, 64 ms 
and 256 ms. In some GRBs with few pulses separated by a long AT the value of T90 which goes from 5% to 95% of 
the counts may not include a small peak before or after the AT yielding a value of AT > T90. This effect is seen in 
(d). Small corrections were applied to the pulse timing locations of the crosses that took into account the truncation 
of the distributions by the limited time resolution. The lines are the best fits to the median values and the power law 
indices are listed in Table 1. The timing values were obtained for pulses that were isolated from neighbouring pulses 
at > 50% whereas all pulses above 5 a were used for AT. 



at 256 ms resolution. The latter GRBs were binned sepa- 
rately despite the overlap with Tg because they are from 
a dimmer sample that was included to extend the range 
in duration to ultra-long GRBs. The GRBs were further 
split into three categories based on the number of pulses N 
> 5cr i.e. 1 < N < 3, 4 < N < 12 and N > 12. The median 
/i and standard deviation a of each lognormal distribution 
was determined ( Quilligan ct al., 2002 ). 



the bursts were combined and the range of the lognormal 
distribution that includes 67% of the values for the rele- 
vant pulse property. The indices of the best fit power-laws 
to the median values in Fig. 1 are listed in Table 1 and all 
have values close to unity. The values of t r , tf , FWHM and 
AT are well correlated with each other dQuilligan et al 



3. Results 
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Table 1. 



Indices of the power law fits to the timing data. 



grams of t r , tf , FWHM and AT versus T90 are pre- 
in Fig. 1 for three categories based on the number 
3S in the GRBs. The crosses in Fig. 1 are not error 
show the range in T90 over which the pulses in 



b it 



Property 



Rise Tin 10 
Fall Tirr e 
FWHM 
Time In erval 



1 < N < 3 



0.98 ±0.07 
0.95 ±0.04 
0.89 ±0.04 
0.97 ±0.06 



4 < N < 12 



0.89 ±0.03 
0.91 ±0.05 
0.87 ±0.03 
0.85 ±0.03 



N > 12 



1.02 ±0.3 
0.72 ±0.1 
0.68 ±0.06 
0.77 ±0.27 
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2002) which explains the similarity of the diagrams when 
plotted versus T90. All four figures are presented to em- 
phasise this result. For the two categories with 1 < N < 3 



by collimated emission with higher average values of T 
close to the rotation axis, where there is lower baryon 



and 4 < N < 12 the power-laws between the median val- 2000). The transparency radius of the fireball r t varies as 



pollution (Rees & Meszaros, 1994; Rees, 1999; Balmonson 



ues of t r , tf , FWHM, AT and T90 extend over a range of 
more than 10 3 and includes GRBs from both sub-classes. 
The GRBs with N > 12 include only those with T 90 > 2 
s because N = 9 is the largest value for GRBs with T90 < 
2 s. The data are displaced from the two catagories of 
GRBs with fewer pulses and display flattening for Tgo < 
50 s which may be caused by a selection effect because of 
pulse amalgamation by the 64 ms resolution of the data. 

The properties of the isolated pulses in GRBs with 
T90 > 2 s, at 64 ms resolution, were compared with the 
preceding (AT - ) and subsequent (AT + ) values of AT. 
The pulse amplitudes (PA) are more anticorrelated with 
AT - than AT + with Spearman rank order correlation co- 
efficients of -0.39 and -0.33. The peak to peak values of AT 
were then corrected for the strong correlation between AT 
and th 



r~2. The profiles of pulses in bursts from off-axis will be 
slowed by (1) the reduced value of T and (2) the longer 
time for faster shocks to catch slower ones before reaching 
rt. There will be fewer shocks to collide outside it and 
hence generate lower luminosity bursts with slower and 
softer pulses because the additional shock amalgamation 
produces a narrower range in T. Indeed this process may 
be a key factor in controlling the spread in T in GRBs. 

The luminosity of B ATSE bursts is not determined but 
several factors imply that GRBs with faster and harder 
pulses are more luminous. The luminosity- variability and 
luminosity-lag correlations for GRBs with known z in- 
fer that the more variable and spectrally harder bursts 
are more luminous (Norris et al., 2000 ; [Salmonson, 2000 



The amended values of the correlation coefficients are 
0.27 and -0.16 with probabilities that the anti-correlations 
arose by chance of 10~ 15 and 10~ 6 . This difference did not 
depend on the isolation level of the pulses which was var- 
ied from 30% to 80%. The FWHM of the pulses were also 
compared with AT~ and AT + and found to be strongly 



Fenimore & Ramirez- Ruiz, 2001 ; [Bchaefer et al., 2001 ; Ioka 



t pulse rise and fall Limes jQifilligan el al., 2002|). \ Nakamura, 2001fc|Schmidt, 2001|). These correlations can 



be explained in an inhomogeneous jet model where the 
more luminous and variable bursts with higher N occur 
close to the axis and the slower, softer and less luminous 



bursts from further off-axis ( |Quilligan et al., 2002| ). In ad- 
dition, there is the strong case for jets in GRBs based 



correlated. The correlations were marginally stronger for 



on th e impressive range of afterglow studies ( Frail et al 



AT" t han AT 1 . 

4. Discussion 

The data in Fig. 1 shows that the median values of t r , 
tf, FWHM and AT depend strongly on T 9 o and also on 
the number of pulses. GRBs with the same durations can 
have a wide variation in the number of pulses and the pulse 
properties. GRBs with a small number of pulses generally 
have slow pulses that are further apart while bursts with 
larger number of pulses have faster and spectrally harder 



200 1| ). The unexpected effect is that the jet has left its 



pulses that are closer together (Norris et al., 2001). In 
both cases the same pattern is observed except that the 
median values of t r , tf, FWHM and AT scale with Tgo. 
Bursts with only one pulse and correlated pulse properties 
would lie on a line of slope unity, while bursts with many 
pulses lie on parallel lines if the pulse properties and AT 
are correlated. The correlated nature of the pulse parame- 
ters and AT provide additional information on the central 
engine and the emission process. 

4.1. Number of pulses, pulse properties and jets 

The separation of GRBs based on N (Fig. 1) suggests a 
kinematic origin because the median values of the pulse 
timing parameters and AT scale in the same way and by 
about the same amount. In the internal shock model, an 
effect of this type may occur in homogeneous jet models 
where the degree of collimatio n and range in T depe nds 
on the mass at the explosion ( Kobayashi et al., 2001 ) or 
in an inhomogeneous jet model by a viewing effect caused 



imprint on the timing profile of the burst. 



4.2. Variation of pulse properties and A T with Tgo 

The power laws in Fig. 1 connect the pulse data in the two 
sub-classes for 1 < N < 3 and 3 < N < 12. The median 
values of the pulse timing parameters and AT scale with 
T90 by about a factor of 10 3 . This result may indicate a 
direct connection between the two sub-classes of bursts 
or that they coincidentally lie on the same power laws in 
the absence of sufficient data to independently determine 
the slopes for the short bursts. The popular progenitors 
of GRBs range from mergers of compact objects such as 
neutron stars and black holes to collapsars and hyper- 
novae in massive stars. A neat feature of most progenitor 
models is that they provide a generic scenario based on 
the formation of a black hole with a massive debris torus 
that is rapidly accreted and energises the jet by neutrino 
transport and magnetohydrodynamic processes. 

The variation of the pulse parameters with Tgo maybe 
caused by the same emission mechanisms and progenitors. 
If GRBs have the same progenitors, the short GRBs must 
have very high values of T to produce the short durations 



(Piran, 1999). The high value of T is close to the upper 
bound allowed by the internal shock model ( [5ari fc Piran 



1997; Lazzati et al., 1999) but well below what can in 



principle be produced in a very low baryon environment 
( Meszaros Sz Rees, 1997 ). An alternative interpretation is 



that there are two types of progenitors with pulses that 
coincidentally lie on the same power laws. The progeni- 
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tors that arise from massive stars may not be capable of 
producing short GRBs. The short bursts may come from 
mergers of neutron stars or neutron stars and black holes 



(Ruffert & Janka, 199E) where the time structure in the 
bursts reflects the interaction of thin relativistic shocks 
with r ~ 100 and duration of ~ 50 ms. In this two progen- 
itor scenario the high values of T can be reduced by about 
an order of magnitude. In addition the central engine and 
the environment before the photosphere may combine to 
smooth the shocks in a way that is a function of Tgo- 



4.3. Correlations between pulses 

There are a number of new results on pulses, from short 
and long GRBs, that constrain the emission process. These 
include (1) the distributions of values of t r , tf, FWHM, 
pulse amplitude (PA) and pulse area of isolated pulses are 
not random but have lognormal distributions (|Quilligan 
|et al., 2002| ; |McBreen et al., 200l| ; |Nakar fe Piran, 200 1| ), 
(2) the values of AT are lognormally distributed with a 
power-law excess of long time intervals or a Pareto-Levy 
tail for GRBs with T 90 > 2 s flQuilligan et al., 2002| ; 



McBreen et al., 2001; Nakar fc Piran, 2001), (3) there is 



an anticorrelation between the values of PA and FWHM 
(Ramirez-Ruiz & Fenimore, 2000; Quilligan et al., 2002| ), 



(4) there is a positive correlation between the values of the 
PA that extend over many pulses (Quilligan et al., 2002; 
McBreen et al., 200l[ ), (5) there is also a positive corre- 
lation between the values of AT flQuilligan et al., 2002 : 
McBreen et al., 2001), (6) the PA of isolated pulses in 



GRBs with Tgo > 2 s and 64 ms resolution were compared 
with AT~ and AT + because of a strong correlation be- 
tween the pulse width and AT~ in the galactic superlumi- 
nal source GRS 1915+105 where outbursts are powered by 



accretion into a black hole (Belloni et al., 1997) and a cor 



relation between long quiescent periods in GRBs and sub- 
seque nt intervals of emission ( Ramirez-Ruiz fc Merloni 
2001 ) . The PAs were found to be anticorrelated with AT - 
and AT + with a slightly larger anticorrelation for AT~. 
Any additional correlations between the pulse properties 
and AT~ would be expected to be much weaker than in 
GRS 1915+105 because a similar intermittent pattern in 
the outbursts from the central engine would have to per- 
sist after shock interactions generated the outbursts. The 
combination of the six results make the pulses in GRBs 
quite unique. In addition long pulses are asymmetric and 
reach their maximum earlier in higher energy bands while 
shorter pulses tend to be more symmetric with negligible 
time lags between energy channels (Norris et al., 1996). 

The three results discussed above seem to favour an 
internal shock model powered by a hyper-accretion pro 



cess into a newly formed black hole ( 


Ruffert & Janka, 


1999; Popham et al., 1999; 


Zhang & Fryer, 2001 


). The 



accretion process has been modelled and is sensitive to 
the rate at which material piles up around the black hole 
and the accretion time scale of the particles. The accre- 
tion rate might provide the overall control on the process 



that eventually generates the correlated outbursts from 
the jet. The pulse amplitudes and time intervals can be 
coupled because higher accretion rates cause larger out- 
bursts that are closer in time, while lower accretion rates 
produce smaller and slower events that are further apart. 
It appears that variation in the rate of accretion, the thick- 
ness of the relativistic shocks and the viewing angle of 
the jet may be key factors in accounting for the observed 
durations and pulse properties in GRBs. There are many 
uncertain factors that influence the process. These include 
the viscosity of the particles, the mass and angular mo- 
mentum of the disk, the mass and spin of the black hole 
and the energy extraction and collimation process. 
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